Dressed-quarks and the nucleon's axial charge 
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The nucleon's axial charge, gA, expresses features that are both fundamental to the strong inter- 
action and crucial to its connection with weak interaction physics. We show that dynamical chiral 
symmetry breaking (DCSB) suppresses the axial-charge of a dressed-quark, g A , at infrared mo- 
menta. Since this effect disappears as chiral symmetry is restored, one may argue that gA vanishes 
with the restoration of chiral symmetry because no nucleon bound-state survives the associated 
transition. The suppression of g q A is shown to be part of an explanation for a 25% reduction of gA 
from its nonrelativistic quark-model value. Critical too, however, is the presence of dressed-quark 
angular momentum within the nucleon. The value of g q A depends on the kernels of the gap and 
Bethe-Salpeter equations. We find that incorporation of essentially nonperturbative effects asso- 
ciated with DCSB into these kernels inflates the value relative to that obtained at leading-order 
in a widely used truncation of QCD's Dyson- Schwinger equations. Such corrections also affect 
the nucleon's axial radius. In both cases, however, agreement with experiment will require similar 
improvements to the Faddeev kernel and associated interaction current. 

PACS numbers: 12. 38. Aw, 14.20.Dh, 12.15.-y, 12.38.Lg 



I. INTRODUCTION 

The prototypical weak interaction is nuclear /3~-decay, 
which explains the instability of neutron-rich nuclei and 
proceeds via the transition 

n -> p + e~ + v e . (1) 

The first attempt at its explanation [l|, Q was based on a 
contact current-current interaction, modulated by a con- 
stant @|: G F = 1.166 x KT 5 GcV~ 2 . Elcctroweak gauge 
theory replaces the contact interaction by exchange of a 
heavy gauge boson and produces the tree-level expres- 
sion (Gp/y/2) = g 2 /(8M^), where M w ~ 80GcV is 
the gauge-boson's mass and g is a universal dimension- 
less coupling; namely, it is the same for all interactions 
between gauge-bosons, leptons and current-quarks. 

Neutron /3-decay and kindred processes play a cru- 
cial role in many domains, e.g.: Big-Bang nucleosyn- 
thesis, constraining the abundance of deuterium; su- 
pernovae explosions, producing a vast amount of en- 
ergy through neutrino production; testing the Standard 
Model, placing constraints on extensions via low-energy 
experiments; and in practical applications, such as 1 C- 
dating and positron emission tomography. Notwith- 
standing its widespread importance, a connection be- 
tween the coupling, g, that describes weak processes in- 
volving current-quarks and that between weak bosons 
and the dressed-quarks produced by nonperturbative in- 
teractions in QCD, the strongly interacting part of the 
Standard Model, has not been elucidated. Attempts to 
do so are described in Refs. |4l4l"lj|. 

The transition in Eq. ([T]) may be studied via the quark- 
based axial- vector matrix element 

A 5m(p/'P*) = (p(p/> A /)I"757^Kp*> A»)) ( 2 ) 



where Pij, Xij are, respectively, initial/final momenta 
and helicities associated with the initial-state neutron 
and final-state proton. If one assumes isospin symmetry, 
then Poincare covariance entails that this matrix element 
is completely described by two form factors 



A Z(Pf'Pi) = u p(Pf' X f) 



«75; 



where q = pf — Pi, Ga{q 2 ) is the nucleon's axial-vector 
form factor, Gp(q 2 ) is its induced pseudoscalar form fac- 
tor and Mn is the average nucleon mass0 The quan- 
tity of primary interest herein is the nucleon's nonsinglct 
axial-charge 

G A {q 2 



9A 



0). 



(4) 



The axial charge also has a relation to nucleon spin 
physics. Given our assumption of isospin symmetry, then 

(p(P/> A /)l «757M d \n(pi,Xi)) 

= 2 (p(Pf, A/) | W757m u ~ d757 M d \p(p h A,)) . (5) 



In the forward scattering limit; i.e., pj 



Pi 



p, with 



A/ = A l = A, A • p = 
frame 

2M N A M (q t ) 
(it) 



~, then in the infinite- momentum 

{pip, A) | 9757^9 |p(Pj A)) , (6) 

f dx [Aq(x) + Aq(x)} , (7) 
Jo 



1 We use a Euclidean metric: {7 M ,7„} = 2<5 M „; 7^ = 7 M ; 75 : 
74717273, tr[757 M 7„7 P 7 (T ] = -4e M „ pCT ; = (i/2)[7 M ,7„]; a-b ■■ 
Yli=i a i^i' and timelike => P 2 < 0. 
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where Aq(x) = qf(x) ~ q±{x) is the light- front helicity 
distribution for a quark q carrying a fraction x of the 
proton's light-front momentum. This measures the dif- 
ference between the light-front number-density of quarks 
with helicity parallel to that of the proton and the den- 
sity of quarks with helicity antiparallel. The connection 
between Eq. ([5]) and helicity will not be surprising, given 
the relationship that may be drawn between the matrix 
structure 757^ and the Pauli-Lubanski four-vector; and 
it follows that 

g A = Z A ( dx [Am(i) + Au(x) - Ad(x) - Ad(x)] ; (8) 



namely, the nonsinglet axial charge measures the differ- 
ence in the light-front frame between the u- and d-quark 
contributions to the proton's helicity 0, [T^]. (Owing 
to the axial-vector Ward-Takahashi identity, the renor- 
malisation constant for the axial-vector vertex satisfies 
Za = Z2, with Z2 discussed below.) 

The induced pseudoscalar form factor, Gp(q 2 ), holds 
its own fascinations, owing particularly to its connection 
with pion-nucleon interactions. Fundamental to the char- 
acter and strength of such interactions is dynamical chiral 
symmetry breaking (DCSB), the phenomenon responsi- 
ble for both 98% of the visible mass in the universe and 
masslessness of the chiral-limit pion [l4|. The existence 
of such a pion entails 



2M 



N 



NN ! 



(9) 



where /° is the pion's leptonic decay constant and g® NN 
is the pion-nucleon coupling constant, where the super- 
script "0" indicates a quantity evaluated in the chiral 
limit. 

Using a Gordon identity, chiral-limit axial-vector cur- 
rent conservation at the nucleon- level (g M A|"(p/,pi) = 0) 
delivers the Goldberger-Treiman relation: 



i\ r° „ a — f° „° 
m n 9a ~ Jtt 5tt 



NN ' 



(10) 



This identity has some curious implications. For exam- 
ple, in the absence of DCSB, /" = and hence no pseu- 
doscalar meson couples to the weak interaction |l5|. It 
then follows from Eq. (JTUJ) that if a nucleon exists with a 
finite, nonzero mass in a universe without DCSB, g A = 
for that nucleon; i.e., such nucleons, too, decouple from 
the weak interaction. (We will subsequently return to 
this.) In these circumstances then g A appears to serve 
as an order parameter for DCSB and a nonzero value of 
OA signals the presence of in-hadron quark condensates 



II. DRESSED-QUARKS 

Our goal is to elucidate a connection between gA and 
the strong physics of dresscd-quarks, which arc described 



in QCD by the gap equation: 

Sj\p)^Z 2 {z 1 -p + mY) 

+Z, jT g 2 D^(p - q )^S f (q)^Tl(q,p), (11) 

where: / denotes the quark's flavour; D^ lL , is the gluon 
propagator; r£, the quark-gluon vertex; J d , a sym- 
bol representing a Poincare invariant rcgularization of 
the four-dimensional integral, with A the regularization 
mass-scale; mj m (A), the current-quark bare mass; and 
Zi j 2(C 2 ) A 2 ), respectively, the vertex and quark wave- 
function renormalisation constants, with £ the renormal- 
isation point. 

The gap equation's solution is the dressed-quark prop- 
agator, 

S f (p) = l/\i 7 -pA f (p 2 ,e) +B f (p 2 ,C 2 )} , (12) 
= Z f (p 2 X 2 )/[^-P + M f (p 2 )}. (13) 

The mass function, M/(p 2 ), is independent of the renor- 
malisation point; and the renormalised current-quark 
mass, 



mj = Z m (C, A) m^ m (A) = Z^Z 2 mf 



(14) 



wherein Z4 is the renormalisation constant associated 
with the Lagrangian's mass-term. The rcnormalisation- 
group invariant current-quark mass may be inferred via 



lim 



1 , P 
- In — 5- 

2 A 2 



QCD 



M f (p 2 ) . 



(15) 



where 7 m = 12/(33 — 2iVj). The chiral limit is 

m/ = 0. (16) 

Chiral symmetry and its breaking pattern in QCD are 
expressed in the following axial-vector Ward-Takahashi 
identity: 



P M rg(fc; P) + i [m/(0 + m g (C)} n s (k; P) 
Sj 1 (k + )i-f 5 + i7 5 S'~ 1 (fc_) , 



,/'.</ 1 



(17) 



I5 are, respectively, amputated axial- 



where and 
vector and pseudoscalar vertices. They connect an out- 
going quark of flavour / and an incoming quark of flavour 
g, with total momentum P = pi + pr and relative mo- 
mentum k = (1 — r))pi + rjpf 1 where r\ G [0, 1], and hence 



k + = pf = k + r)P, k- 



Pi 



k — (1 — rf)P. Owing 



to Poincare covariance, no observable can legitimately 
depend on 77; i.e., the definition of the relative momen- 
tum. N.B. Equation (JTTJ) is modified for flavourless pseu- 
doscalar mesons and this leads to important differences 
in their behaviour, which are discussed in Ref. (22l |. 

The vertices relevant to /3~-decay are r^J, T^ d but 
with our assumption of isospin symmetry we can ignore 
the flavour labels and consider the diagonal elements 
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~pu— t 
1 5^ 



r 5/1! rr 



The axial-vector vertex then 



has the general form [23j 

r 5AI (fc; P) = 7 5 hMk; P) + ■ kG R (k; P) 
-a^KH R (k;P)]+t 5fi {k;P) 



P, 



p2 



■2/ ff r,(fc;P) 



(18) 



where: Fr, Gr, H r and r5 M (fc;P) are regular in 



the neighbourhood of [P 2 



0; P M r 5M (fc;P) 



P a PpM a p(k;P), with M a p(k;P) a matrix-valued func- 
tion; and the pion's Bcthc-Salpctcr amplitude is 

T w (k;P)=-f 5 [iE„(k;P)+-fPF„(k;P) 

+ fc • P 7 • fc G w (fc; P) + a^k^H^k; P)] . (19) 

Combining now Eqs. (jTTJ) , (|12l) . (|T71) - ((T9^) and working in 
the chiral limit, one may derive |23| the following quark- 
level Goldbcrgcr-Trciman relations: 



f°E^(k;0) = B°(k 2 ), 
F R (k;0)+2f°F^k- : 0)^A°(k 2 ), 

G R (k; 0) + 2f°G r (k; 0) = ^A\k 2 ) 



H R (k;0)+2f%H w (k;0) = 0. 



(20) 
(21) 

(22) 

(23) 



These identities arc of critical importance in QCD. The 
first exposes the fascinating consequence that the solu- 
tion of the two-body pseudoscalar bound-state problem 
is almost completely known once the one-body problem is 
solved for the dressed-quark propagator: the relative mo- 
mentum within the bound-state is identified unambigu- 
ously with the momentum of the dressed-quark. This last 
fact emphasises that Goldstone's theorem has a pointwisc 
expression in QCD. It is difficult to overestimate its im- 
portance for Standard Model physics. 

The remaining three identities show that a pseu- 
doscalar meson Goldstone boson must contain compo- 
nents of pseudovector origin. Some of the important 
corollaries of this result are exposed in Refs. [24| - |28| . 
Herein, however, we reveal additional novel consequences 
of Eqs. 



III. AXIAL CHARGE OF A DRESSED-QUARK 

Consider the dresscd-quark-axial-vector vertex, 
Eg. (TTgj) . Only F R (k;P), the function associated with 
the Dirac structure 757^, possesses an ultraviolet 
divergence in QCD perturbation theory. Notably, in 
Landau gauge the renormalised amplitude F R = 1, 
up to next-to-leading-order perturbative corrections: 
one-loop corrections vanish. (This may be derived fol- 
lowing Ref . [29| . ) All other functions in the axial- vector 
vertex are power-law suppressed in the ultraviolet. In 
perturbation theory, therefore, the quantity 



expresses the distribution of a current-quark's axial- 
charge. It remains perturbatively close to unity. (The 
impact of other components in Eq. (|18j) is canvassed in 
Sec. El They do not materially affect our discussion.) 

Nonperturbatively, however, the situation is very 
different, as may readily be illustrated. To this 
end, consider the symmetry-preserving regularisation of 
a vector x vector contact-interaction detailed and em- 
ployed in Refs. Hii| As elucidated 
therein, in rainbow-ladder truncation^ such an interac- 
tion produces results for low-momentum-transfer observ- 
ables that are practically indistinguishable from those 
generated by more sophisticated interactions, such as 
that explained in Refs. [H|, [35|. The consequences of 
Eq. (|2"Tj) are dramatic in this context. With the single 
parameter determining the interaction strength chosen 
small, am/ir < 0.4, then DCSB is absent and 



q DCSB 



(25) 



where "CN" denotes contact interaction. On the other 
hand, with arR/7r ~ 1; namely, chosen commensu- 
rate with contemporary estimates of the zero-momentum 
value of a running-coupling in QCD (HI, I36T4381 ] . one 
obtains A°(k 2 ) = 1, M°(k 2 ) = M° = 0.358 GeV, 
M°F®(k; 0) = 0.46, all fc-independent with a contact in- 
teraction, and f® = 0.1 GeV, so that 



P° (k; 0) = 1 -2/0 F°(fc;0) = 0.74. 



(26) 



Thus the quantity associated with the current-quark's 
axial-charge is markedly suppressed in the infrared owing 
to the nonperturbative phenomenon of DCSB. 

To allay any concern that this outcome might be 
model specific, we compared Eq. (|26[) with the value pro- 
duced by the most sophisticated rainbow-ladder interac- 
tion available [Hj], which is detailed in App. lA ll In this 
case one naturally finds a fc 2 -dependent form for g q A and 
obtains 9a rl (^ 2 = 0) = 0.81 at a realistic value for the 
light-quark current-mass. 

We are also able to compare these results with that 
produced by the most complete kernels available for the 
gap- and Bethe-Salpeter equations [H|. These kernels, 
described briefly in App. lA 21 and denoted subsequently 
by "DB," incorporate essentially nonperturbative effects 
associated with DCSB, such as a dressed-quark anoma- 
lous magnetic moment [40l - |42^ . and yield 



9at 



,(0)= 0.87= 1.06 g ARL (0) 



(27) 



The infrared suppression is thus seen to be a generic fea- 
ture of the axial- vector vertex. Its impact is far-reaching 
since it will influence, e.g.: the leptonic radiative decays 



g q A (k 2 ) :=F R (k;P = 0) 



(24) 



2 Rainbow-ladder is the leading-order in a systematic and 
symmetry-preserving truncation scheme for QCD's Dyson- 
Schwinger equations [32I, |33|1 ■ 
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of charged light pscudoscalar mesons; and the nucleon's 
axial charge, as we shall subsequently see. 

We stress that the infrared suppression of the axial- 
vector vertex contrasts markedly with the effect of dress- 
ing on the leading covariant, j^, in the vector vertex. 
In this case the associated scalar function is bounded 
below by unity at (k = 0; P = 0) owing to the vector 
Ward-Takahashi identity. Indeed, with a momentum- 
dependent interaction this scalar function is always en- 
hanced, as illustrated in Fig. 2 of Ref. [39| . 

At this point it is worth emphasising that Poincare 
covariance demands that the general form for a pseu- 
doscalar meson Bethe-Salpeter amplitude possess four 
components; namely, those appearing in Eg. (|19p . In- 
spection of the Bethe-Salpeter equation for pseudoscalar 
mesons shows that a nonzero value for E„ is the force 
behind ^ 0, with the coupling fixed by the DCSB 
mass-scale, which is provided by the dressed-quark mass- 
function, M. To be clear, M ^ in the chiral limit entails 
E n 7^ 0, and together these results require F^ ^ 0. Read- 
ily apparent in the rainbow-ladder truncation, this is a 
general result, independent of the function chosen to rep- 
resent the drcsscd-gluon and the Ansatz for the dressed- 
quark-gluon vertex in the gap equation, Eq. For 
further confirmation, compare the results reported above 
with those in, e.g., Refs. (^IMSfil- 

This explains why the appearance of F n =^ is a nec- 
essary consequence of DCSB. Given an interaction with 
nontrivial momentum dependence, then G n and H v are 
also necessarily nonzero for the same reason. Plainly, 
a complete expression of DCSB is not achieved merely 
by producing nonzero values for the in-pion condensate 
and pion leptonic decay constant. The full structure of 
the Goldstone mode must also be described. Finally, 
positivity of /° guarantees that of ^(fcjO), and hence 
the second term on the left-hand-side of Eq. (|2"Tj) is pos- 
itive. This means that Fji(k 2 ; P = 0) is bounded above 
by A°(k 2 ) and approaches this function from below as 
k 2 -> oo. (This is illustrated in Fig. 8 of Ref. Q.) It 
does not, however guarantee Fr(0; 0) < 1. That is a con- 
sequence of the dynamics which produces the Goldstone 
pion and sets the mass-scale for DCSB. 

It is now a propos to reconsider the role of qa in connec- 
tion with DCSB. In chiral-limit Dyson-Schwinger equa- 
tion (DSE) studies, chiral symmetry restoration and de- 
confinement are coincident no matter which control pa- 
rameter is varied0 This supports a view that DCSB and 
confinement are intimately related; and we expect that in 
the presence of some agent which undermines the inter- 
action strength required for DCSB, confinement is also 
lost. Under these conditions /° = and consequently 



F R (k 2 ;0) = A°(k 2 ), following from EqflHJ. Should such 
circumstances correspond to a domain whereupon none 
of the interactions in the Standard Model is strong, then 
both functions will be unity up to perturbative correc- 
tions. On the other hand, suppose that strong corre- 
lations remain after chiral symmetry restoration, such 
as may be in a putative strongly-coupled quark-gluon 
plasma, then Fn(k 2 ; 0) = A°(k 2 ) > 1: i.e.. both functions 
are actually enhanced above unity [4(| |47|. Evidently 
then a connection between the restoration of chiral sym- 
metry and g A vanishing is not driven by changes at the 
level of the axial- vector dressed-quark vertex. 

Consider now that the identity in Eq. (|10jl holds so 
long as chiral symmetry is dynamically broken and a nu- 
cleon exists with nonzero and finite mass, even under 
conditions that place the theory in the neighbourhood 
of /° = + . Given that g q A > 1 in these circumstances, 
a vanishing of g A would require extraordinary and pre- 
cise cancellations amongst the terms that constitute the 
nucleon's axial-charge matrix element; i.e., between the 
various contributions arising from the angular momen- 
tum correlations within the nucleon's Faddeev ampli- 
tude. Owing to the power of symmetries in quantum 
field theory, this is conceivable but nevertheless improb- 
able. Given the preceding discussion it is more likely 
that the chiral-limit relationship g A — > is connected 
with dissolution of the nucleon bound-state at a point 
of coincident chiral symmetry restoration and deconfine- 
ment. A realisation of this phenomenon is illustrated for 
the scalar and pseudoscalar meson sector in Sec. IV of 
Ref. 0| . The conjecture may be tested using modern 
Faddeev equation treatments of the nucleon. 

A vanishing of gA entails that the right-hand-side of 
Eq. (jHJ is zero. This expression is normally described as 
expressing the difference in the light-front frame between 
the u- and d-quark contributions to the proton's helicity. 
How can that vanish? One is here considering the chiral 
limit. Absent a DCSB mechanism, a chiral limit theory 
with massless quarks separates into two distinct, non- 
communicating theories: one for positive helicity states 
and another for negative helicity. Each sub-theory has 
identical interactions and hence each will produce the 
same quark number distributions, labelled, however, by 
opposite helicities. Since there is no mechanism in the 
total theory that can flip helicity, the number of positive 
helicity states will always match the number with nega- 
tive helicity. Hence the result <?a = is achieved because 
each of the four terms in Eq. jS]) vanishes individually, 
irrespective of whether or not they are associated with a 
bound-state. 



3 See, e.g., the discussions in Refs. [4(j|47j |. for which it is important 
to note that light-quark confinement is not connected in any 
known manner with the static-quark potential. It can instead be 
related to marked differences between the analy tic p roperties of 
coloured and colour-singlct Schwingcr functions |48|| . 



IV. QUARK MODELS AND SA 

Related to constituent-quark model phenomenology, 
Eqs. ([26]) . (|27 |) are curious. It is textbook knowledge (see, 
e.g., Ref. |50|) that constituent-quark models with spin- 
flavour wave- functions based on SU(6) symmetry pro- 
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duce the following axial-charge of the nuclcon: 



Q 
9 A =2,9a 



d 3 x 



u 2 (x) - \v 2 (x) 



l ~\J d 3 xv 2 (x) 



(28) 
(29) 



where g A is the axial-charge of a constituent-quark, and 
u(x), v(x) are, respectively, the upper and lower compo- 
nents of the nucleon's constituent-quark wave-function. 
Plainly, in a nonrelativistic model, v(x) = and gA = 
(5/3) g®, so that reproducing the empirical value of gA = 
1.27 requires g^ NR — 0.76. This value compares well 
with those in Eqs. (|2"6")l , (f2"7) . Of course, the origin of the 
empirical value of gA is more complicated but nonpertur- 
bative dressing of g q A plays a part. 

A full explanation is suggested by Eq. (f2"9"| , which has 
two key features. As we have described above, the first is 
dressing of the axial- vector vertex, an effect that modifies 
the strength with which a dressed-quark couples to the 
W-boson. 

The other is indicated by the second term within the 
parentheses in Eq. (f2"9")l : 



Cv = - I d 3 xv 2 (x). 



(30) 



This represents the appearance of P-wave quark orbital 
angular momentum in a relativistic constituent-quark 
model. 

In a quantum field theory such as QCD, the nucleon 
is properly described by a Poincare covariant Faddeev 
equation |5l|. In this context, may be reinterpreted as 
signifying the impact of correlations within the nucleon's 
Faddeev wave-function, which possesses S-, P- and D- 
wave dressed-quark orbital angular momentum compo- 
nents in the nucleon's rest frame. In the presence of 
DCSB, such correlations are strong. For example, the 
5-wave-only contribution to the nucleon's normalisation 
is just 60% [10, HI1 and it is known that altering the 
strength of quark orbital angular momentum correlations 
within the nuclcon can materially affect gA [53| . Within 
this framework, therefore, the empirical value of gA em- 
bodies the outcome of interference between dressing the 
quark-I-F-boson vertex and angular momentum correla- 
tions within the nucleon's Faddeev amplitude. This is 
not too surprising given the connection between gA and 
the u- and d-quark hclicity distributions, expressed in 
Eq. Notably, however, the magnitude of the sup- 
pression of g q A and the strength of orbital angular mo- 
mentum correlations in bound-state wave functions are 
both driven by DCSB. 



4 The canonical normalisation constant for the nucleon's Faddeev 
amplitude is equivalent to requiring that the nucleon's Dirac form 
factor is unity at zero momentum transfer. 



The latter is readily seen from Sec. Ill of Ref. [54[. In 
the absence of DCSB, the amplitudes F n , G^, in 
Eq. (fT9|) vanish identically in the chiral limit, as does 
M(p 2 ) in Eq. (|13p . It follows that any correlation that 
survives is described by a Bethe-Salpeter wave-function: 
XTx{k\P) = S(p)T 7r (k;P)S(p) cx 75, a purely 5-wave 
structure in rest-frame kinematics. 

These causal relationships emphasise again, following 
identical conclusions drawn from computations of the 
pion and nucleon valence-quark distributions pil [30L [55l — 
|58I |. that understanding parton distribution functions 
(PDFs), as opposed to merely parametrising them, rests 
upon grasping the nature of DCSB in QCD. It exposes 
the potential gains to be made in hadron physics by shift- 
ing theoretical focus from modelling PDFs to their well- 
constrained computation. In this connection it should 
be borne in mind that the first three non-trivial PDF 
moments, which is the maximum that can be obtained 
from numerical simulations of lattice QCD (57| . are insuf- 
ficient for PDF reconstruction: more than ten moments 
are required in order to constrain the large-x exponent 
to better than 10% [H. 



V. FADDEEV EQUATION AND g A 

It is known that Faddeev equation models can be con- 
structed to reproduce the empirical value of gA l60l | , 
unifying it in the process with other nucleon observables. 
However, such studies employed axial- vector vertices that 
do not respect Eqs. (J2H)— (E3J) - This is mended in Ref. H|, 
which solves all elements of the problem - the gap, Bcthc- 
Salpetcr and Faddeev equations - in rainbow-ladder trun- 
cation. That study, however, produces gA RL — 0.99(2), 
underestimating the empirical value by 22%. Natu- 
rally, with M% and /° near to their experimental values, 
Eq. (fTU|) entails that Ref. Q underestimates g^N n by a 
similar amount. 

The magnitude of the error is typical of rainbow-ladder 
truncation in those channels for which it is known and 
understood a priori to be adequate. In the sector of 
light-quark vector and flavour nonsinglct pseudoscalar 
mesons, over an illustrative basket of thirty-one calcu- 
lated quantities tabulated in Ref. [6l| , the truncation de- 
livers a standard-deviation of 15% in the relative error 
between experiment and theory Q. 

Part of the remedy to this quantitative error lies in go- 
ing beyond the leading-order truncation when solving the 
gap and Bethe-Salpeter equations. This is now possible 
in a symmetry-preserving manner [62| . as indicated by 
Eq. (|27|) . Indeed, we have solved for the dressed-quark- 
axial- vector vertex using the kernels described in Ref. [39| 
and recapitulated in App. lA 2[ which are essentially non- 
perturbative, incorporating effects of DCSB that were 
not previously possible to express. These kernels clarify 
a causal connection between DCSB and the splitting be- 
tween vector and axial-vector mesons, and expose a key 
role played by the anomalous chromomagnctic moment of 
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FIG. 1. Selected functions in the axial- vector vertex, Eq. (|31[1 : 
rainbow-ladder (RL) result cf. that obtained with DCSB- 
improved kernels for gap and Bethe-Salpeter equations (DB). 
Curves, all dimensionless: solid FP B and dashed 

F RL. 

very- 
long-dashed 2M§ B Ff B and dot-dashed 2Mf L F 3 HL ; and dot- 
dash-dash (-4Mf b2 F? b ) and dot-dot-dash (-4M| l 2 F 7 rl ). 

(mE b 



0.36 GeV and M§ L = 0.41. 



dressed-quarks [40h42L l63l |64J in determining the values 
of observable quantities. 

The general form of the transverse part of the axial- 
vector vertex, all that contributes directly to g\, is 

r^(A; P) = 75 [-tfFi - i 7/ h • Pk ■ PF 2 + T^a up k p F 3 
+ [fc^ 7 • P + i^o vp k v P p \F± - ifcj-fc ■ PF 5 



+k p L -/ ■ Pk ■ PF e + fc^ 7 ■ kF 7 
+k^a vp k v P p k ■ PF 8 ] , 



(31) 



where: {Fi\i = 1,...,8} are scalar functions of 
(k 2 ,k ■ P,P 2 ) that are even under k ■ P — s- (— k ■ P); 



P 2 



1; T„ 



[5 pv - P p P v /P 2 ], T pv + L pv = 8 pv ] and 
a^ L = T pu a v . Since Eq. (|3 1 [) may simply be obtained from 
Eq. psp through contraction with T^, we have the fol- 
lowing correspondences: F\ -H> Fr, F 7 f-> Gr, F 3 <-> if^. 

Given that a dressed-quark anomalous chromomag- 
netic moment produces a large dressed-quark anomalous 
magnetic moment [42[ , one should at least expect that F$ , 
with its similar tensor structure, is significantly altered 
when proceeding beyond rainbow-ladder truncation. In 
fact, all the scalar functions are materially modified on 
a domain < \k\JME ;$ 5, where Me is the Euclidean 
constituent-quark mass, {Me} = {\/s\s = M 2 (s),s > 
0}: -Fi.2,3,5,6,8 magnitudes are enhanced, with F5 also 
changing sign; and -F47 magnitudes are suppressed. In 
Fig.[T]we illustrate the response of each of those functions 
appearing in Eqs. (f2T|) — ([25 ]> . 

With at least eight quantities reacting markedly to im- 
provements in the DSE kernels, it is natural to seek a sin- 
gle measure that can illustrate the plausible consequences 
for gA- To this end we consider 



FIG. 2. With the elements defined in association with 
Eq. (J33J), the ratio \g q A DB (k 2 )\/\g q A RL (k 2 )\ {solid curve). The 
straight line with dotted outliers represents the band 1.064 ± 
0.003. 



where, at each value of p 2 > 0, the Euclidean spinors 
satisfy 7 -pu{p) = <; p u(p), u(p)j-p = u(p)? p , <; p = M(p 2 ). 
Focusing on the case (fc • P = 0, P 2 = 0) and using the 
appropriate Euclidean- Gordon identities, Eq. (f32|) yields 
an axial-charge distribution, which is complex for the on- 
Euclidcan-mass-shell dressed-quarks: 

gf{k 2 ) = FrC^O) + F 3 (, 2 ;Q) . (33) 

With this kinematic arrangement, no other functions 
from Eq. (f3"Tj) contribute. We reiterate and emphasise 
that Eq. (|33|) is not an observable but rather an illustra- 
tive artifice: a simple quantitative measure of the im- 
pact of terms in Eq. (|3ip on the infrared behaviour of the 
axial-charge of a dressed quark. 



In Fig.[Uwe plot the ratio \g% qDB (k 2 )\/\g^ nl " (k 2 )\, 
which is the single measure we sought. It assumes the 
value 1.064 ± 0.003, consistent with Eq. ([2"7|). If one sup- 
poses that corrections to the Faddeev equation result for 
gA, arising from bettering the rainbow- ladder computa- 
tion of the dressed-quark axial-vector vertex, can simply 
be estimated by rescaling the axial-charge of a dressed- 
quark; viz., g A RL 9a DB = l-06g^ RL , then one infers a 
value of = 1.05(2) from Ref. Q. Consistent with stud- 
ies of the nucleons' electromagnetic form factors [EH, this 
is an important but modest improvement. 

The modest size of the improvement is good be- 
cause the utility of rainbow-ladder truncation would have 
been much reduced if the magnification were too large. 
Notwithstanding this, the constant rescaling probably 
underestimates the effect, given the structure apparent 
in Fig.[TJ One would better reckon the correction by 
building an Ansatz for T§ p , consistent with the algebraic 
constraints and numerical results we have elucidated, and 
employing that in a Faddeev equation computation of the 
nucleon's weak and strong form factors. One should bear 
in mind, however, that correcting the gap and Bcthe- 
Salpctcr equation kernels is not the complete picture. 



E„RL 1 



u( Pf )r^(k-p)u( Pl ) 



(32) 
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FIG. 3. P dependence of the functions appearing in Eq. ((33} . 
Curves, all dimensionless: solid F[' B {k 2 = 0;P 2 ) and dashed 
F? L (k 2 = 0;P 2 ); and long-dashed 2M% B F£ B {k 2 = 0; P 2 ) 



and dot-dashed 2MjpF 3 iiL 
and Mj? L = 0.41.) 



(* a 



0:P" 



0.36 GeV 



The Faddeev equation kernel and associated interaction 
current should also be modified. These modifications, 
too, must affect qa- 

In association with this we note that our dressed ker- 
nels do not contain pieces that might reasonably be de- 
scribed as corresponding to meson-cloud effects. Consid- 
ered analysis of such contributions is expected to further 
increase the value of qa by < 10% [H, \6]§. They can be 
added to our kernels and interaction current without con- 
cern for overcounting and hence their effect may also be 
explored. The role played by a meson-cloud in forming qa 
is, in fact, much discussed. A contemporary effective field 
theory perspective may be traced from Ref . [H, [68[ ; that 
within lattice-QCD from Refs. [^4TTj : and that within 
models of nucleon structure from Refs. [1^, [7(| . 

Improvement of the gap and Bethe-Salpeter kernels 
should also affect the result for the nucleon's axial radius, 
ta, quite simply because the rainbow-ladder truncation 
is unable to explain the location of the ai-meson pole in 
the axial-vector vertex, whereas the DCSB-corrected ker- 
nels, App. lA2[ resolve this longstanding problem (39|. In 
underestimating the mass of the a\ meson, the rainbow- 
ladder truncation overestimates the contribution to r a 
from the associated pole. It is likely, therefore, to overes- 
timate this radius or, equally, understate the mass-scale, 
ttia, that characterises evolution of the nucleon's axial 
form factor in the neighbourhood of P 2 = 0. 

In Fig. [3] we depict the P 2 dependence of the functions 
in Eq. (|33f . Once again, the effects of improving the ker- 
nels are measurable. An indication of its impact is the 
ratio of mass-scales that characterise monopole fits to 
the Fi functions in the figure: a DB /a RL = 1.65. This 
matches well with the expectation just described. It is 
curious, however, because Ref. || reports to^ l = 1.28(6), 
which is already at the upper limit of values inferred from 
experiment (reviewed in Ref. 71|): ttia € [1.0,1.3]. This 



corresponds to a value of ta at the lower limit of ex- 
periment. It is relevant here to note that corrections to 
the Faddeev kernel and associated interaction current can 
plausibly magnify correlations within the nucleon and 
their impact on interactions, as the discussion of Fig. [I] 
shows they do for quark-antiquark systems. Such effects 
would serve to increase ta- 



VI. EPILOGUE 

The nucleon's axial charge, gA, expresses features that 
are both fundamental to the strong interaction and cru- 
cial to its connection with weak interaction physics. It is 
thus important to understand how its strength originates 
within QCD and is thereby connected with dynamical 
chiral symmetry breaking (DCSB), the source of more 
than 98% of visible mass in the Universe. 

In pursuing this goal we demonstrated that DCSB sup- 
presses the axial-charge distribution of a dressed-quark, 
g A , at momenta k < Me, where Me ~ 0.4 GeV is 
the mass-scale associated with DCSB. Conversely, quark- 
level Goldberger-Treiman relations indicate that g\ ~ 1 
in the absence of DCSB. This result, combined with the 
nucleon's Goldberger-Treiman relation, led us to a view 
that gA vanishes with the restoration of chiral symmetry 
because no nucleon bound-state survives the associated 
transition. 

Consistent with inferences based on constituent-quark 
models, we found a suppression of g\ to be part of an 
explanation for the value of the nucleon's axial charge. 
Critical, too, however, is the presence of dressed-quark 
angular momentum correlations in the nucleon's rest- 
frame Faddeev wave- function, and hence in (almost) ev- 
ery frame as a result of Poincare covariance. (It would be 
an exceptional Poincare transformation that transferred 
an observer to a frame in which every vestige of orbital 
angular momentum was eliminated.) 

The Poincare covariant Faddeev equation is a natu- 
ral means by which to describe the structure of the nu- 
cleon bound-state. As one of QCD's Dyson-Schwingcr 
equations (DSEs), it is a critical element in a unified 
symmetry-preserving explanation of meson and baryon 
properties. At leading-order in the most commonly used 
truncation (namely, rainbow-ladder approximation), this 
approach produces a value of gA that is 22% smaller than 
experiment. We explained that this is a good result given 
the established level of accuracy that one may expect at 
leading-order 

Complementing this, we argued that well-constrained 
improvements to the kernels of the gap and Bcthe- 
Salpeter equations, which incorporate essentially nonper- 
turbative corrections to the rainbow-ladder truncation, 
increase g q A and are therefore likely to improve the DSE 
result for gA- Such corrections also affect the nucleon's 
axial radius. In both cases we saw that agreement with 
experiment will require similar improvements to the Fad- 
deev kernel and interaction current. 
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In closing it is worth reiterating that one should gener- 
ally expect a ~ 15% mismatch between experiment and 
results obtained in the internally consistent application of 
rainbow-ladder truncation in those channels for which the 
truncation is most reliable; namely, ground-state light- 
quark vector and flavour nonsinglet pseudoscalar mesons, 
and the nucleon and A ground-states. Uniformly pre- 
cise agreement would indicate serious deficiencies in the 
method - a misuse of degrees-of-freedom, for example - 
and diminish materially its capacity to provide insights 
into strong QCD. 

The truncation's simplicity is a strength. One can tol- 
erate such modest disagreement with experiment when 
the result is ready computation of a diverse array of phe- 
nomena, their semiquantitative connection with funda- 
mental elements in QCD, and enabling of a clear identi- 
fication of familial relationships between them. 

Notwithstanding these features, improvement is now 
possible and necessary. Fuller incorporation of DCSB 
into bound-state kernels and interaction currents will en- 
able the better informed feedback between experiment 
and theory that is necessary to understanding: confine- 
ment and DCSB in QCD; the nature of connections be- 
tween them; and how they affect observables ranging 
from elastic and transition form factors to parton dis- 
tribution functions. 



form of Qm,{k 2 ) determines whether confinement and/or 
DCSB are realised in solutions of the gap equation. 
The interaction in Ref. 03 is 



Sir 2 



De~ s l u 



87T 2 7m 5(S) 



ln[r+(l + S /A| CD ) 



(A2) 



where: 7m = 12/(33-2%), N f = 4, A QCD = 0.234 GcV; 
t = c 2 — 1; and J (s) = {1 — exp(— s/[4m 2 ])}/s, m t = 
0.5 GcV. With Duj = constant, light-quark observables 
are independent of the value of a; € [0.4, 0.6] GeV. We 
used Duj = (0.8 GeV) 3 and uj = 0.5 GeV. 

In the rainbow-ladder truncation, in the isospin sym- 
metric limit, the inhomogeneous axial-vector Bethc- 
Salpctcr equation is 

g(k - q) D^(k - q) laX5f i(q; P)lp , (A3) 

dq 

where X 5» = S(q+)T£°(q; P)S(q-). 

Regarding renormalisation of the gap and inhomoge- 
neous Bethe-Salpeter equations, we follow precisely the 
procedures of Refs. [43|, uM an d use the same renormali- 
sation point; i.e., £ = 19 GcV. A current-quark mass of 
ra c = 3.4 MeV produces = 0.136 GeV. 
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Appendix A: Interaction kernels 
1. Rainbow- ladder 

In a sophisticated rainbow-ladder DSE study the 
model input is expressed in a statement about the nature 
of the gap equation's kernel at infrared momenta, since 
the behaviour at momenta k 2 > 2 GeV 2 is fixed by per- 
turbation theory and the renormalisation group [43l . l72j . 
In Eq. (jll[) . this amounts to writing (k = p — q) 

Z l9 2 D^(k)T v {q,v) = k 2 Q(k 2 )D^ c (k) lv 
= [k 2 Q m {k 2 ) + 47rcip Q cD(fc 2 )] D*?(k)-y v , (Al) 

wherein D l ™ c (k) is the Landau-gauge free-gauge-boson 
propagator; a P QCD(fc 2 ) is a bounded, monotonically- 
decreasing regular continuation of the perturbativc-QCD 
running coupling to all values of spacclike-fc 2 ; and 
Qm(k 2 ) is an Ansatz for the interaction at infrared mo- 
menta: Qi K (k 2 ) < <5 pQC d(£: 2 ) Vfc 2 > 2 GeV 2 . The 



2. DCSB-improved kernel 

The DCSB-improved kernel is specified by a dressed- 
quark-gluon vertex and a Bethe-Salpeter kernel deter- 
mined therefrom. 

In the gap equation, Eq. (jTTj) . we use [3!| 

Z l9 2 D pa {t)T a (q,q + t) = g(t 2 )D i ^(t)Z 2 t a (q 1 q + t), 

(A4) 

with Q from Eq. fA2]), w = 0.5 GeV but Duo = 
(0.52 GeV) 3 , a change required to ensure the dressed- 
kernels produce physical observables which match those 
obtained in rainbow-ladder truncation when that is reli- 
able; 



^(pi.pa) - r^(p 1)P2 ) + rr>i,p 2 ) ; 



(A5) 



iTl c { PuP2 ) = iZ A (pl,pi) 7 „ + 2^ [i 7 • t AaIpIpD 

+A B (pl,p 2 2 )] , (A6) 
where the first term was introduced in Ref. (74|, with 

xm,pI) = [m) + m)\/% ^(pIpD = mpD- 

0(pi)]/[Pi — pIL 2^ = Pi + P2\ and the anomalous chro- 
momagnetic moment piece is [42| 

rjr(pi,pa) = r^ cm4 (pi,p 2 ) + r^( PllP2 ) , (A7) 

with (k =pi -p 2 , = Sft,, - k^kv/k 2 , := T^ v a v ) 

r^ cm4 = [^1 ■ k + i-yJa vp l v k p )n( Pl ,p2) , (A8) 
<y^k v T 5 {pi,p 2 ) , (A9) 

2r 5 (pi,p 2 ) 



pacms 



M{plp 2 ) ' 



(A10) 
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r 5 = r)A B (pi,pl), V = 0.65 [2S| ; and M(x,y) = [a 
M(xf +y + M(y) 2 ]/(2[M(x) + M(y)]). 

The inhomogeneous Bethe-Salpcter equation is 



A" 



z 2 g(k~ q )D { ;™(k~q)— la s( q+ ) 



dq 

xT Bli { q] P)S{q-)^p(Q-,k-) 

f \ a 

+Zi g 2 D al) {k- q) — la S f {q + ) 

Jdq Z 



' dq 
\a 

x—A 5flP (k,q; P), 



(All) 



where the four-point function As^p is completely de- 



fined [62j via the quark self-energy and hence the 
quark-gluon vertex, T^. Crucially, A 5m/ j satisfies a 
Ward-Takahashi identity whose solution provides a 
symmetry- preserving Ansatz consistent with T ll . We use 

2A 5/3(m) = [fp(q+,k+) +7 5 f /3 (g_,fc_)7 5 ] 

x s-Hk + ) + s-i(-k^ {k;P) 

+ r 5M (q;P) irT 



S-l(-q + ) + S-Hq-) 
x [ 75 f p (q+, k + ) l5 + f>(g_, (A12) 

Regarding renormalisation, here, too, we follow the 
procedures of Refs. [H, [73| and use the same renormali- 
sation point; i.e., ( = 19 GeV. A current-quark mass of 
m c = 3.7 MeV produces = 0.138 GeV. 
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